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Abstract 

We treat high-energy neutrino production in GRBs. Detailed calculations of photomeson neu- 
trino production are presented for the collapsar model, where internal nonthermal synchrotron 
radiation is the primary target photon field, and the supranova model, where external pulsar-wind 
synchrotron radiation provides important additional target photons. Detection of ^ lOTeV neu- 
trinos from GRBs with Doppler factors > 200, inferred from 7-ray observations, would support the 
supranova model. Detection of < lOTeV neutrinos is possible for neutrinos formed from nuclear 
production. Only the most powerful bursts at fluence levels > 3 x 10 -4 erg cm -2 offer a realistic 
prospect for detection of v^. 

PACS numbers: 14.60.Lm, 95.85 Pw, 95.85.Ry, 98.70.Rz, 98.70Sa 
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Two leading scenarios for the nature of the sources that power lone-duration GRBs are 
the co llapSar j and s_ (SA) Q m „de ls . The co r e of a maS sive star collapses direct* 
to a black hole in the collapsar model, but only after an episode of neutron-star activity 
in the SA model (^J. This delay means that a shell of enriched material, which could 
lelp explain observations of X-ray features, surrounds the GRB source in the SA model 
3]. The colliding shells thought to operate in the collapsar model take place at distances 
T 2 cAt ~ 3 x 10 15 r3 00 At(s) cm from the central ejection source, where T = 300r 30 o is a 
typical wind Lorentz factor and At is the time between shell ejection events Q. GRB pulses 
of ~ 0.1-10 s durations and separations are typical [5], but even for At ~ 1 ms, the shell 
collisions take place far outside the photospheric radii of likely GRB stellar progenitors. The 
most important radiation field for photomeson neutrino production in the collapsar scenario 
is thus the internal synchrotron radiation 

The presence of a pulsar wind (PW) within an expanding supernova remnant (SNR) shell 
in the SA model means that the radiation environments in the collapsar and SA models are 
vastly different |3<j. To maintain stability against prompt collapse to a black hole, the neutron 
star must rotate with periods near 1 ms |2j. The neutron star radiates electromagnetic, 
leptonic, and possibly also hadronic energy in the form of a powerful relativistic 'cold' 
magnetized wind consisting of quasi-monoenergetic e + -e _ pairs and ions. The ordered flow 
of the wind is disrupted at the wind shock formed within the SNR shell. Here we consider 
only the emission from quasi-monoenergetic leptons that are injected at the PW shock 
which provide, primarily through synchrotron losses, the main external photon target for 
photomeson interactions. The quasi-thermal radiation field of the SNR shell is neglected, 
as well as the nonthermal synchrotron and Compton-scattered radiations from particles 
accelerated at the PW shock. 

We find that the external lepton-wind synchrotron radiation alone improves prospects 
for neutrino detection by orders of magnitude over values calculated in a standard external- 
shock/GRB blastwave scenario []| that explains well the phenomenology of GRBs with 
smooth fast-rise/ slow decay 7-ray light curves. The presence of the external field can increase 
the number of detectable neutrinos by an order of magnitude or more over a colliding shell 
scenario that is generally invoked to explain highly variable GRB light curves. 

Here we report detailed calculations of photomeson neutrino production for the collapsar 
and SA models that take into account nonthermal proton injection followed by photomeson 
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energy loss, which is computed numerically by adapting our photo-hadronic model for blazar 
jets j8|. We also calculate the corresponding cutoff spectral energy e 77 due to 77 pair- 
production attenuation. High-energy neutrino observations of very bright GRBs, especially 
if combined with data from the Gamma-ray Large Area Space Telescope (GLAST) or other 
gamma-ray detectors, can test the collapsar and SA models, as we now show. 

The detection efficiency in water or ice of ultrarelativistic upward-going muon neutrinos 
with energies e„ = 10 14 ei4 eV is P Vil = lCr 4 e^ 4 , where x = 1 for e 14 < 1, and x — 0.5 
for ei4 > 1 j9(. For a neutrino ffuence spectrum parameterized by v<§ v = lO -4 0-4 e i4 erg 
cm -2 , the number of detec 
A v = 10 10 A 10 cm 2 is therefore 

V ^v n A r,r 0-4^10 



cm 2 , the number of detected with a km-scale v detector such as IceCube with area 
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e 14 , tor eu > 1 . 

(1) 

For a v^ v spectrum with a v ~ 0, the number of to be expected are N u ~ 1.20_ 4 ^4 1O (1 + 
| In e^) for ei4 < 1, and N v ~ 1.20_ 4 ylio/ A /ii4 for ei 4 > 1. If the nonthermal proton energy 
injected in the proper frame is comparable to the radiated energy required to form GRBs 
with hard X-ray/soft 7-ray ffuences > 10~ 4 ergs cm -2 , then extremely bright GRBs are 
required to leave any prospect for detecting „„ with km-scale nentrino detectors. About ,2-5 
GRBs per year are expected with hard X-ray/soft 7-ray ffuence > 3 x f 4 ergs cm 2 [10]. 

Because GRB blast waves are believed to accelerate nonthermal electrons, it is probable 
that they also accelerate nonthermal hadrons. Coincidence between the GRB power radiated 
at hard X-ray/soft 7-ray energies within the GZK radius and the power required to account 
for super-GZK 10 20 eV) cosmic rays suggests that GRBs might be the progenitor sources 
of ultra-high energy (> 10 19 eV) cosmic ravs, cosmic rays above the ankle of the cosmic ray 
spectrum, and GeV-TeV cosmic rays jll| \\. 

In our calculations, we inject protons with a number spectrum oc e~ 2 at comoving proton 
energies e p > 300r 30 o GeV up to a maximum proton energy determined by the condition 
that the particle Larmor radius is smaller than both the size scale of the emitting region and 

n 

the photomeson energy- loss length [8J. The observed synchrotron spectral flux in the prompt 
phase of the burst is parameterized by the expression F(u) oc v^ 1 (v / Vbr) a , where hv^ T = 300 
keV, a = —0.5 above z/& r , and a = 0.5 when 30keV < hv < hv^, T . At lower energies, a = 4/3. 
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The observed total hard X-ray/soft 7-ray photon fluence <&tot — tdur Jo°° duF(u), where td ur 
is the characteristic duration of the GRB. In our calculations we assume a source at redshift 
z = 1, and let $ tot = 3 x 10~ 5 erg cm -2 . Two or three GRBs should occur each month above 
this fluence level. 

We inject a total amount of accelerated proton energy E' = 47id 2 L § t ot8~ 3 (l + z)~ l into 
the comoving frame of the GRB blast wave. Here 5 is the Doppler factor and di is the 
luminosity distance. The energy deposited into each of N sp light-curve pulses (or spikes) is 
therefore E' = E' /N sp ergs. We assume that all the energy E' is injected in the first half 
of the time interval of the pulse, which effectively corresponds to a characteristic variability 
time scale t var = t rfur ./2A^ sp . The proper width of the radiating region forming the pulse 
is AR' = t var c5/(l + z), from which the energy density of the synchrotron radiation can 
be determined ||S|. We set the GRB prompt duration td ur = 100 s, and let N sp = 50, 
corresponding to t var = 1 s. The magnetic field is determined by assuming equipartition 
between the energy densities of the magnetic field and the electron energy. 

Fig. 1 shows the optical depth r 77 for 77 pair production attenuation as a function of 
observer photon energy, calculated for 5 = 100,200, and 300. Curves 1, 2, and 3 show the 
77 opacity from internal synchrotron radiation only, corresponding to the collapsar scenario. 
The rapid decrease of r 77 with increasing 5 is explained by the rapid decrease of the energy 
density of the internal radiation in the comoving frame as u' syn ~ L' syn /27iR' 2 c oc {vF v ) ohs / 5 & , 
which implies a rapid decline (approximately, r 77 cx 5~ 4 , depending precisely on the radiation 
spectrum) of the r 77 opacity in the synchrotron field. Relativistic flows with 5 ^ 100 
are required to explain observations of > 100 MeV 7 rays with EGRET \l2\. GLAST 
observations may imply more stringent limits on 5. 

An external radiation field given by the expression vL u oc v 1 ! 2 exp(— z//z/ exi ), with hu ext = 
0.1 keV, is assumed to be present in the SA model [3|. The intensity of this field is determined 
by the assumption that the integral power L ext = J °° L v du is equal to the power of the 
pulsar wind L pw w (10 53 erg)/td e iay, assuming that a total of ~ 10 53 erg of pulsar rotation 
energy is radiated during the time tdeiay (which is here set equal to 0.1 yr) from the rotating 
supramassive neutron star before it collapses to a black hole. The energy hv ext ~ 0.1 keV is 
the characteristic energy of synchrotron radiation emitted by electrons (of the pulsar wind) 
with Lorentz factors 7 pMJ ~ 3 x 10 4 in a randomly ordered magnetic field of strength w 10 G. 
The radius R is determined by assuming v = 0.05c is the mean speed of the SNR shell, and 
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the external photon energy density oc Lj 2nR 2 . The thin solid curve is the combined opacity 
in a SA-model pulse when 5 = 100. In this case, the pulses will be highly attenuated above 
« 300 MeV. 

Fig. 2 shows the total fluences expected from a model GRB with N sp = 50 pulses. The 
thin curves show collapsar model results at 5 = 100, 200, and 300. The expected numbers of 
Up that a km-scale detector such as IceCube would detect are N v = 3.2 x 10~ 3 , 1.5 x 10~ 4 , 
and 1.9 x 10~ 5 , respectively. (The effect of neutrino flavor oscillations could reduce these 
numbers by a factor rs 2.) There is no prospect to detect from GRBs at these levels. 
The heavy solid and dashed curves in Fig. 2 give the SA model predictions of N v = 0.009 
for both 5 = 100 and 8 = 300. The equipartition magnetic fields are 1.9 kG and 0.25 kG, 
respectively. The external radiation field in the SA model makes the neutrino detection rate 
insensitive to the value of 5 (as well as to t var > 0.1 s, as verified by calculations), but there 
is still little hope that a km 3 detector could detect such GRBs. 

Neutrino production efficiency would improve in the collapsar and SA models if t var ~ 
1 ms and N sp = 5 x 10 4 to provide the same total fluence. We obtain N v = 0.027, 0.012, and 
0.0046 for 5 = 100, 200, and 300, respectively. Such narrow spikes are, however, then nearly 
opaque to gamma rays, with r 77 = 1 at E 1 = 80 MeV when 5 = 200, and at = 500 
MeV when 5 = 300. When t var ^> 1 s, the GRB blast wave is optically thin to its internal 
synchrotron radiation at multi-GeV energies. In the absence of an external radiation field, 
the produced fluxes of neutrinos are not however detectable. 

Neutrino detection from GRBs can be assured only if the number of background counts 

i- f poo patm / Q \ 2 -ex, 

B{> e u ) ^ J dQJ dtj de l P Ufl A u ~ 5xlO~ 8 iMio f — ) J dxx ~^ +X ~ 2 < L ( 2 ) 

Here the chosen time window has duration t w = 10 2 t2 s, 9 is the half-opening angle, and the 
cosmic-ray induced atmospheric neutrino background flux at the nadir is F^ tm w 8 x 10~ n e7/ 
ergs cm" 2 s -1 sr -1 , with (3 = 1.7 for e i4 < 1, and (3 = 2 for eu > 1 [9]. This expression places 
a lower bound on the energy above which background counts may be neglected, namely that 
the neutrino energy e u > e min = 3.6(t 2 Ao) a59 (#/l°) L18 GeV in order to have B < 1. Note 
that IceCube's angular resolution is predicted to be ~ 1° at TeV energies, 0.5° at 50-100 TeV 
energies, and 0.6° at PeV energies jy]. Therefore a detection of only 2 neutrinos during the 
time interval of the prompt phase of a GRB, especially if at energies e v > 10 TeV predicted 
by photohadronic interactions, will be highly significant. 



The importance of low-energy (10 GeV - 10 TeV) sensitivity in neutrino experiments is 
connected with the possibility to probe jet models with nuclear interactions. One possible 
scenario is a beam-on-target model where the GRB protons pass through and occasio nally 
collide with particles in a dense target, namely the SNR shell in the SA scenario [13] . 
Neutrinos created from the decay of mesons formed in interactions with SNR shell particles 
with mean atomic weight A are beamed into an angle no smaller than the opening angle of 
the GRB outflow. Thus the neutrino fluence is at most comparable to the hard X-ray/soft 
7-ray fluence of a GRB multiplied by the conversion efficiency r\ p = t sh /3t p A of protons to 
neutrinos, again assuming as before that the GRB hadronic energy is determined by $>tot- 
The factor of 3 accounts for the fraction of secondary energy emitted in the form of v^. Here 
the available time for the interaction of relativistic protons passing through the SNR shell 
is t s h = 3.3 x 10 4 /_ii?i6 s, assuming a uniform shell with a characteristic thickness 0.1/_i-R 
and radius R = 10 16 i?i6- The pA — > 7T ±,CI energy-loss time scale t p A in the stationary frame is 
given by t~\ = K p n A a p AC ~ 4.3 x IQ~ 7 m S N r/ {A 1 ^ R\ & f -i) s _1 , given an inelasticity K p ~ 0.5 
and a strong interaction cross section <j p a — 30A 2 / 3 mb. Thus r] p ~ 5 x lO~ 3 msNRRii A^ 1 / 3 . 
These neutrinos are formed at energies ~ 0.05Tm p c 2 ~ 15r3oo GeV. Neutrinos formed 
through beam/target interactions are therefore relatively low energy, and are not expected 
from GRBs unless Ri$ <C 1, in which case the SNR shell would be very Thomson thick. 

Another nuclear interaction scenario is where relativistic particles in the GRB blast 
wave collide with other particles in the ejecta. The relativistic particles are formed ei- 
ther by sweeping up and isotropizing particles captured from the external medium 
or by accelerating particles at external or internal shocks. The z/^-production efficiency 
Vpp — t'ava/Stpp, where t' ava < R/Tc is the available time in the comoving frame, t'^ 1 = 
n' p a pp c = E iso cr pp c/ (4TcR 2 ARTm p c 2 ) is the inverse of the proper-frame nuclear-collision time 
scale, and E iso = 10 52 E 52 ergs is the apparent isotropic GRB energy release. Because 
AR' ~ rct var , we find r\ pp < 6 x 10" 4 -E , 52/[-Ri6r| 00 /f: 1 , ar (s)]. Secondary neutrinos are formed 
at energies ~ 0.05T 2 m p c 2 = 4T 2 00 TeV. Only GRBs with large values of E 52 or small values 
of r 30 o, Rw, and t var (s) provide reasonable efficiencies for production, and these would ex- 
hibit strong 77 attenuation with comparable fluences in and reprocessed electromagnetic 
radiation. 

Thus detection of a few neutrinos by a km-scale detector generally requires an extremely 
bright event at the level reaching ~ 10~ 3 erg cm -2 . Note in this regard that a model recently 
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proposed for neutrino emission from the pp interaction of PW protons with the SNR 
shell during the time td e i ay between supernova and GRB events neglects the expected flux 
of accompanying electromagnetic radiation, such that a pre-GRB source in the SA scenario 
would be an extremely bright source significantly before the GRB event, which would be 
hard to miss in all-sky X-ray and 7-ray observations. A recent study of neutrino production 



during the prompt GRB phase in the SA model 



171 ] finds similar estimates for the number 



of expected neutrinos which, however, are predicted to be mostly at much higher enegies, 
> 1000 TeV. Moreover, photomeson interactions of protons on the internal synchrotron 
radiation field as well as 77 attenuation are not studied in Ref. [3, 3|, which is important 
to discriminate between different GRB models. 

Another way to improve high-energy neutrino and UHECR model production efficiency 
is to assume that the radiating particles in the GRB blast waves are dominated by hadrons. 
As Fig. 2 shows, the integrated level of fluence and therefore of secondary gamma-rays for 
the assumed injection of relativistic protons is at the level ^ 3 x 10 -6 erg cm -2 ~ 0.1$^. 
It would therefore still be possible to assume acceleration of protons with a power up to 
a factor « 10 higher than in Fig. 2, which would increase the number of expected for 
IceCube to ~ 0.1 in the SA model. The detection of a few neutrinos would then be possible 
in the prompt phase of GRBs at the flux level &tot ~ 3 x 10~ 4 erg/cm 2 . It is important that 
from photomeson interactions only neutrinos above 10-20 TeV are to be expected, while a 
detection of lower energy neutrinos would significantly favor models with nuclear interactions 
in either the SA or the collapsar scenarios. 

In this letter, we have shown that the collapsar model is less efficient for neutrino pro- 
duction than the supranova model, as it lacks the external PW synchrotron radiation field. 
Neutrino predictions with the collapsar model are very sensitive to 5, and are most favor- 
able when the ~ 100 MeV - GeV 77 opacity is large. An inverse relation between the 
number of detectable neutrinos and cutoff spectral energy is found for a collapsar model. 
Neutrino production efficiency is relatively insensitive to 5 in the SA model, though it does 
depend on model parameters of the PW synchrotron emission. There is no prospect for 
neutrino detection in the collapsar model when 5 200. More optimistic estimates from the 
viewpoint of detecting GRB neutrinos could be found in proton-dominated GRB models. 
Without this hypothesis only the brightest GRBs can be expected to be detected with both 
high-energy 7-ray and neutrino detectors. Detection of high-energy neutrinos from GRBs 
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would therefore have far-reaching impact on the collapsar and SA scenarios for GRBs, the 
hypothesis that GRBs sources can be powerful accelerators of ultrarelativistic CRs, and for 
our understanding of the origin of GRB radiation. 
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FIG. 1: Opacity r 77 to pair-production attenuation for 7 rays which are detected with observer 
energy given on the abscissa. The radiating region moves with Doppler factor 5 with respect to the 
observer. The internal synchrotron radiation field is typical of GRBs observed with BATSE, and 
the injected energy corresponds to a measured fluence of 3 x 10 -5 ergs cm -2 , which is radiated in 50 
equal pulses. The heavy and light dashed curves give the 77 opacity through the external radiation 
field over the size scale of the supernova remnant and the pulse emitting region, respectively. The 
thin solid curve gives the total opacity for a model pulse in the external shock/SA model. 
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FIG. 2: Energy fluence of photomeson muon neutrinos for a model GRB. The thin curves show 
collapsar model results where only the internal synchrotron radiation field provides a source of 
target photons. The thick curves show the 5 = 100 and 300 results for the SA-model calculation, 
which includes the effects of an external pulsar wind radiation field. 
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